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ABSTRACT 

The cross sections for the v—d and v—d reactions are calculated for the incident energy up 
to E^ = 170 MeV, with the use of a phenomenological Lagrangian approach. We assess and 



O ' improve the reliability of the employed calculational method by examining the dependence 

^ ■ of the results on various input and approximations that go into the calculation. The main 



points of improvements over the existing work are: (1) use of the "modern" NN potentials; 
(2) use of the more accurate nucleon weak-interaction form factors; (3) monitoring the 



2 ■ strength of a vertex that governs the exchange-current contribution, with the use of data on 

the related process, n + p ^ d + 'j. In addition to the total cross sections, we present various 
differential cross sections that are expected to be useful for the SNO and other experiments. 
In the low energy regime relevant to the solar neutrinos, the newly calculated total cross 
sections essentially agree with the existing literature values. The origins of slight differences 
found for higher energies are discussed. The ratio between the neutral-current and charged- 
current reaction cross sections is found to be extremely stable against any variations in the 
input of our calculation. 



I. INTRODUCTION 

The neutrino-deuterium reaction^ have been studied extensively over the past decades 
||lHl^- Recent detailed studies are strongly motivated by the proposal and successful start 
of the Sudbury Neutrino Observatory (SNO) [jl5|, which uses a large underground heavy- 
water Cerenkov counter. One of the primary goals of SNO is to study the solar neutrinos 
by monitoring three reactions occurring in heavy water: (i) v — e scattering, z/g + e^ — *■ 
z/e + e~; (ii) the charged-current (CC) reaction, v^ ^ d -^ e~ + p + p; (iii) the neutral- 
current (NC) reaction, u^ + d ^ u^ + n + p, where z/^, stands for a neutrino of any flavor. 
The unique feature of SNO is its ability to register the CC and NC reactions separately 
but simultaneously. Since the NC reaction measures the total flux of the solar neutrinos 
(regardless of their flavors), SNO experiments offer valuable information about the nature 
of possible neutrino oscillation. SNO is also capable of monitoring astrophysical neutrinos 
the energy of which extends well beyond the solar neutrinos energy regime, a prominent 
example being supernova neutrinos. Obviously, in interpreting experimental results to be 
obtained at SNO, the accurate knowledge of the zz — c? reaction cross sections is a prerequisite. 
Although the z/ — e scattering cross section is readily available from the standard model, 
estimation of the neutrino-deuteron reaction cross sections requires a detailed examination 
of the structure of two-nucleon systems and their responses to electroweak probes. 

In describing the current theoretical situation regarding the u—d cross sections, it is useful 
to consider the u — d reactions in a broader context of the general responses of two-nucleon 
systems to electroweak probes. A highly successful method for describing these responses 
is to consider one-body impulse approximation terms and two-body exchange-current terms 
acting on non-relativistic nuclear wave functions, with the exchange currents derived from 
a one-boson exchange model. In a modern realization of this approach |jl6|-0 , the vertices 



^ When convenient, we use the word "neutrino" and the symbol 'V in a generic sense, referring 
to both neutrinos and anti-neutrinos. 



characterizing relevant Feynman diagrams are determined, as much as possible, with the use 
of the low-energy theorems and current algebra. Some coupling constants are inferred from 
models (the quark model, SU(3), SU(6), etc.). In the present work we refer to this type of 
formalism as the phenomenological Lagrangian approach (PhLA). This formalism has been 
used extensively for electromagnetic processes in two nucleon systems |]I^|-pT||. The reported 



good agreement between theory and experiment gives a strong hint of the basic soundness 
of PhLA. This method has also been applied to two-nucleon weak-interaction processes such 
as muon capture on the deuteron Il8| j2^j23[| , the pp-iusion reaction [^,|^], and the u — d 
reactions. For the muon capture, the calculated capture rate agrees reasonably well with 
the experimental value, again rendering support for the basic legitimacy of PhLA. (For the 
pp fusion there is unfortunately no data available.) 

For the neutrino-deuterium reactions, the most detailed study within the framework of 
the impulse approximation (lA) has been done by Ying, Henley and Haxton (YHH) ||10|| , 
while the most elaborate PhLA calculations including the exchange- current effects as well 
as the lA terms have been carried out in IISHnil, and the latest status is described in [112 



to be referred to as KN.|^ In the solar neutrino energy regime, the cross sections given in 
KN are slightly larger than those of YHH. This difference, however, is mostly due to the 
absence of the exchange-current contributions in YHH. As far as comparison with data is 
concerned, Tatara et a/.'s estimate [§] of o"(z/e + d —^ e~ + p + p) averaged over the Michel 
spectrum of z/g agrees with the result of a stopped-pion-beam experiment [^ within large 
experimental errors (30%). Furthermore, the result of a Bugey reactor neutrino experiment 
p6| agrees, within 10% experimental errors, with the values of cr{iJe + d —^ e'^ + n + n) and 
a{ue + d ^ Ue + p + n) given in KN. Thus PhLA seems to provide a reasonably reliable 
framework for calculating the neutrino-deuteron cross sections. 



^ Ref. [^ also gives a rather detailed account of the relation between these latest calculations 
and the earlier work. 



Meanwhile, a new approach based on effective field theory (EFT) has been scoring great 



success in describing low-energy electroweak processes in the two-nucleon systems p7|-|32 
In particular, the rate of thermal neutron radiative capture on the proton {n + p -^ d + '-/) 
has been calculated in chiral perturbation theory (xPT) and the result is found to be in 



perfect agreement with the data |^^. Butler and Chen [^ and Butler, Chen and Kong |T1 
have recently made extremely elaborate studies oi u — d cross sections for the solar neutrino 
energies with the use of EFT. The results of their EFT calculation agree with those of PhLA 
in the following sense. In an EFT approach, one starts with a general effective Lagrangian, 
Ceff, that contains all possible terms compatible with given symmetries and a given order 
of expansion; the coefficient of each term in Cgff is called the low-energy coefficient (LEC). 
Now, it often happens that some LEC's cannot be fixed by the symmetry requirements alone 



and hence need to be treated as parameters to be determined empirically. In [|T3|,|T4[, the 
coefficient Lia of a four-nucleon axial-current counter term enters as an unknown parameter, 
although dimensional arguments suggest — 6fm^ < Lia < -|-6fm'^. According to [|13|, the 
v — d cross sections obtained in EFT agree with those of the PhLA calculation (YHH or 
KN), provided Lia is adjusted appropriately. The optimal value of Lia is Lia = 6.3 fm^ for 
YHH, and Lia = 1.0 fm^ for KN, reasonable values as compared with the above-mentioned 
dimensional estimates. The fact that an ab initio calculation (modulo one free parameter) 
based on EFT is consistent with the results of PhLA provides further evidence for the basic 
reliability of PhLA. 

Bahcall, Krastev and Smirnov [Q have recently studied in great detail the consequences 
of measurements of various observables at SNO. As input for their analysis, the Ve~d reaction 
cross sections of YHH and KN are used, and the difference between these two calculations 
are assumed to represent la theoretical errors. According to |Q , uncertainties in the v — d 
cross sections represent the largest ambiguity in most physics conclusions obtainable from the 
SNO observables, a feature that again points to the importance of reducing the uncertainty 
in the v — d reaction cross sections. 

In the present article we carry out, within the framework of PhLA, a detailed study of the 



cross sections for the charged- current (CC) and neutral-current (NC) reactions of neutrinos 
and anti-neutrinos with the deuteron: 

i^e + d ^ e^ +P + P, (1) 

ly^ + d^u^+p + n (x = e, /i, or r), (2) 

Ue + d -^ e^ + n + n, (3) 

Px + d^Ux+p + n (x = e, /i, or r) . (4) 

It is our view that, in calculating the low-energy v — d cross sections, EFT and PhLA 
play complementary roles. EFT, being a general framework, is capable of giving model- 
independent results, 'provided all the LEC's in an effective Lagrangian C^ff are predeter- 
mined. At present, however, C^ff does contain an unknown LEG, /^ia-I Meanwhile, al- 
though PhLA is a model approach, its basic idea and the parameters contained in it have 
been tested using many observables. Thus, insofar as one accepts the validity of these tests, 
PhLA has predictive power. It is reassuring that, as mentioned, there is highly quantitative 
correspondence |T^,n] between the low-energy u — d cross sections obtained in PhLA and 



those of EFT within a reasonable range for Lia. In this article we wish to investigate several 
key aspects of PhLA in more depth than hitherto reported. 

Beyond the solar neutrino energy regime, PhLA is at present the only available formalism 



for evaluating the u — d cross sections. The EFT calculation in ||13| , ^4| , by design, "integrates 
out" all the degrees of freedom but that of the heavy baryon. The nature of this so- 
called "nucleon-only" EFT limits its applicability to very low incident neutrino energies 
(typically the solar neutrino energies) .0 On the other hand, there is no obvious conceptual 



In principle, however, it is possible to fix Lia using a parity-violating electron scattering exper- 



iment. Hjlj 

One can hope to extend the applicability of EFT to higher energies by including the pion degree 
of freedom explicitly via xPT. An ab initio calculation based on xPT for the i' — d reactions is yet 



obstacle in using PhLA in an energy regime significantly higher than that of the solar 
neutrinos. Therefore, once the reliability of PhLA is tested at low energies by comparison 
with experimental data or with the results of EFT, it is rather natural to use PhLA for 
higher energies as well. In this sense, too, EFT and PhLA seem to play complementary 
roles (at least, in the current status of the matter). 

Our main goals here is to assess and improve the reliability of the PhLA calculation of the 
u — d reaction cross sections, by carefully examining the dependence of the results on various 
input and approximations that go into calculations. The main points of improvements in 
this work over the existing estimates are: (1) Use of the "modern" NN potentials; (2) Use 
of the more accurate nucleon weak-interaction form factors; (3) Monitoring the strength of 
the irNA vertex that governs by far the dominant exchange-current contribution, with the 
use of data on the related process, n + p ^ d + '-/. A second practical goal of this paper 
is to provide detailed information about the various differential cross sections for the v — d 
reactions. Although the total cross sections are well documented in the literature, there 
have not been systematic descriptions of the differential cross sections. We therefore discuss 
in detail the energy spectrum, angular distribution and double differential cross sections of 
the final lepton in the CC reaction and also the energy spectrum and angular distribution of 
the final neutron in the NC reaction. It is hoped that, the detailed information given here 
on these differential cross sections will be useful in analyzing SNO and other experiments. 

In the low energy regime relevant to the solar neutrinos, our results are found to be in 
essential agreement with those of KN. Based on these and additional results described in this 
article, we shall deduce the best estimates of theoretical errors in the v — d cross sections. 
For higher energies, the present calculation gives v — d total cross sections larger than those 
of KN by up to 6%; we shall discuss the origin of this variance. 

The organization of the rest of this article is as follows. After giving in Sec. II a brief 



to be done. 



account of the general framework of our PhLA, we describe in Sec. Ill the calculational 
details, including the multipole expansion of the nuclear currents, and expressions for the 
cross sections for v — d reactions. The numerical results are presented in Sec. IV, and 
discussion and summary are given in Sec. V. Some kinematical formulae necessary for 
calculating phase space integrals are given in Appendix. 

II. FORMALISM 

We are concerned with the v /v — d reactions listed in Eq.(P-Eq.(^. The four-momenta 
of the participating particles are labeled as 

v/v{k) + d{P) -. l{k') + Nr{p\) + N2{p'^), (5) 

where £ corresponds to e^ for the CC reactions [Eqs.([^),(^], and to z/ or z/ for the NC 
reactions [Eqs.(0),(^]. The energy-momentum conservation reads k + P = k' + P' with 
P' = p[ +P2) cind we denote a momentum transfer from lepton to nucleus by q'^ = k'^ — k''^ = 
P'^ — P^. In the laboratory system to be used throughout this work, we write 

F = {E,, k), k'^ = {E'„ k'), P^ = (M,, 0), P'^ = (P'°, P'), (1^ = (^, Q)- (6) 

The interaction Hamiltonian for semileptonic weak processes is given by the product of 
the hadron current {J\) and the lepton current (L'^) as| 



CC _ Gf cos Oc 
for the CC process and 



fdx[J^^{x)L\x) + h. c] (7) 



rjNC ^F 



V2 



Jdx [J^^\x)L\x)+h. c] (8) 



^ Throughout we use the Bjorken and Drell convention for the metric and Dirac matrices, except 
that we adopt the Dirac spinor normahzed as n^n = 1. 



for the NC process. Here Gp = 1.166 x 10~^ GeV~^ is the Fermi couphng constant, and 
cos^c = 0.9749 is the Cabibbo angle. 
The lepton current is given by 

L\x)=i,iixh\l-^')Mx), (9) 

and its matrix element is written as 

l^ =< k'\L^{0)\k > =ui{k')-f^{l--f^)u^{k) for //-reaction, 

= Vp{k)'y^{l — 7^)f j-(A;') for //—reaction . (10) 

The hadronic charged current has the form 

J^^{x) = V,H^) + A^{x), (11) 

where Vx and Ax denote the vector and axial-vector currents, respectively. The superscript 
+ (— ) denotes the isospin raising (lowering) operator for the z/(z/)-reaction. Meanwhile, 
according to the standard model, the hadronic neutral current is given by 

Jf ^(a;) = (1 - 2 sin^ 9w)V^ + AI-2 sin^ OwV^, (12) 

where 9w is the Weinberg angle with sin^ 9w = 0.2312. Vx is the isoscalar part of the vector 
current, and the superscript '3' denotes the third component of the isovector current. In the 
present case the hadron current consists of one-nucleon impulse approximation (lA) terms 
and two-body meson exchange current (MEX) terms. Their explicit forms are described in 
the next subsections. 

A. Impulse approximation current 

The lA current is determined by the single-nucleon matrix elements of Jx- The nucleon 
matrix elements of the currents are written as 

<N{p') I V^{0) I N{p) > = u{p')[fy^x + ^^^Apg1r±w(p), (13) 

<N{p') I A±(0) I N{p) > = u{p')[fAlxl' + fpl'qx]rMp) , (14) 



where Mat is the average of the masses of the final two nucleons. For the third component 



of the isovector current, we simply replace r"^ with y- For the isoscalar current 

<Nip') I 1^,^(0) I N{p)> =u{p')[fy^, + z^a,,q^]\u{p). 
The non-relativistic forms of the lA currents are given by 

i 



yfAfi^x 



VfA\ 



X 



■^lAfiV^ 



Ma{x 



V/aA^ 



VU^ 



, '2Mn 



2Mn 2Mn 

o-i ■ {p't + Pi) - 
, fp 



2M, 



N 



2Mn 
V{V-(Ti)]Tt6{x-ri 



^ /y - 6{x - Vi 



Uf^ 



Pi + Pi , fv + l 



+ 



A/' 



1 



-V X crA - 6(x — Ti). 
2Mn 2 ^ ' 



(15) 

(16) 
(17) 
(18) 
(19) 
(20) 

(21) 



2M^ 

It is useful to rewrite Pi + p[ = q + P ± 2pjv, where the +(— ) sign corresponds to i = 1 
{i = 2), and the derivative operator pj^ should act on the deuteron wave function; in the 
laboratory system we are working in, we have P = 0. 

As for the g^ dependence of the form factors we use the results of the latest analyses in 



3l,|35|: 



Mql) = GD{ql){l+ fipv)il + v)~\ 

fhlii) = GD{q^^){fip - /in - 1 - finV){'^+V) 



-1 



fAiql 



fpiQl) 



-1.254 GAiql), 



2M 



N 



ml~ql 



fAiql), 



fldiql) = GD{ql){f^p + /i„ - 1 + /i„r7)(i + 7]) \ 



with 



GdU 



Qt 



GAiql] 



0.71GeVV 

2 \ -2 

ql \ 



1.14GeV' 



(22) 
(23) 
(24) 
(25) 
(26) 

(27) 
(28) 



where fXp = 2.793, /i„ = —1.913, r] = — ^ffr and m^ is the pion mass. 
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N 



B. Exchange currents 

As mentioned, we use a phenomenological Lagrangian approach (PhLA) to estimate the 
contributions of meson-exchange (MEX) currents. In a PhLA due to Ivanov and Truhhk 
1 17 1, the MEX operators are derived in a hard pion approach |36], in which one exphcitly 



constructs a phenomenological Lagrangian consistent with current algebra, PCAC, and the 
vector meson dominance. This Lagrangian was used by Tatara et al. in their calculations 
ioT fi — d capture and the u — d reactions |P . Meanwhile, the studies by Doi et al. PJ23[ 



indicate that only a small subset of the possible diagrams gives essentially the same results 
as the full set. Based on this experience, we consider here the following types of exchange 
currents. 

Axial-vector current 

The axial vector exchange current, ^mex? consists of a pion-pole term and a non-pole 
part, A'^j^x- Using the PCAC hypothesis, we can express A'^^x i^ terms of the non-pole 
part alone: 

^MEX = "^MEX 2 2 (^ ' ^MEX — ^AmeX,o)- (29) 

We therefore need only consider the non-pole part. For the time component it is known that 
one-pion exchange diagram gives the most important contribution, called the KDR current 
|37|f| The explicit form of the KDR current, with a vertex form factor supplemented, reads 

+(1 « 2), (30) 



As discussed extensively in [38,39], corrections to the KDR current can arise from heavy-meson 
exchange diagrams. We however do not consider those corrections here, since the contribution of 
the KDR current in the present case turns out to be small (see below). 

10 



with r = ri — r2 and iu.,^ = Jq'"^ + rn^. For the space component, we take account of 
the isobar current, A J, that arises from one-pion and one-p-meson exchange diagrams. Its 
exphcit form is 

^^^^{Cq'rt^ + rfi(^i X q')\T^ x T2](^)}(^2 ■ q') 



X 






KiQ 



UJ^. 



\cpq' X (0-2 X q')Tf^ + dpai x (q' x (erg x q'))[Ti x T2]^^^} 



p 
+ (1 ^ 2), (31) 



with Up = \ q'"^ + ^Wp and nip is the mass of p-meson. For the third component of the 
isovector current, we just replace t^" and [ti x t-2\^^^ with ^ and ^"^^^^^^ — , respectively. 
(This prescription will be applied to other exchange currents as well.) The numerical values 
of the pion coupling constants can be determined from low-energy pion-nucleon scattering 
pO| , while the p-meson coupling constants are deduced from the quark model: 

^ = 0.08, coml = 0.188, diml = -0.044, Cpml = 36.2, dp = -\cp. 

4:71 

Furthermore, we assume that the g^-dependence of the vertex form factors, KmN{q^) and 
KmAiq^) {m = 7T,p), are given by K^N{q^) = K^A{q^) = K^{.q^) = ^^^.KpN^q^) = 
Kp/^{q^) = Kp{q^) = ^^, with cutoff masses, A^ = 1.18 GeV and A^ = 1.45 GeV 
Pl| . We use the above-listed values of couphng constants and form factors as our standard 
parameters. 

Vector current 

Regarding the vector exchange currents, we first note that the exchange currents for the 
time component must be small, since the exchange currents for charge vanish in the static 
limit. As for the space component, we take into account pair, pionic, and isobar currents. 
If we adopt the one-pion exchange model for the pair and pionic current and the one-pion 
and one-p-meson exchange model for the isobar current, their explicit forms are given as 

11 



+(1 « 2) , (32) 

g-iq'i-(ri-a;) g-iq'2-(3:-r2) 

X 2 2 ('^i ■ 9'i)('^2 ■ q'2){q'i + q'2), (33) 

Vi(-) = -^tii^V X At (34) 



with uj^i = Jml + q, 



C. Nucleon-nucleon potential 



In PhLA, the nuclear transition matrix elements are obtained by sandwiching the one- 
body lA and two-body MEX currents between the initial and final nuclear wave functions 
which obey the Schrodinger equation that involves a phenomenological nucleon-nucleon po- 
tential. The earlier work PJ^ indicates that, so long as we use a realistic NN potential that 
reproduces with sufficient accuracy the scattering phase shifts and the deuteron properties, 
the numerical results for the v — d cross sections are not too sensitive to particular choices 
of NN potentials. It seems worthwhile to further check this stability for the modern poten- 
tials that were not available at the time of the work described in [§,0. As representatives 
of the "state-of-the-art" NN potentials, we consider in this work the following three: the 
Argonne-tilS potential (ANLV18) [^, the Reid93 potential [^, and the Nijmegen II po- 



tential (NIJ II) [^]. For the sake of definiteness, however, we treat ANLV18 as a primary 
representative. We shall compare our results with those obtained with the use of the more 
traditional potentials. 

D. Monitoring the reliability of the model 

Although, as mentioned, there is by now a rather long list of experimental and theoretical 
work that points to the basic robustness of PhLA calculations, it is desirable to monitor the 

12 



reliability of our model by simultaneously studying reactions that are closely related to the 
V — d reactions and for which experimental data is available. It turns out that the vriVA 
vertex that features in the dominant exchange current for the v — d reaction appears also in 
the np -^ '-yd reaction, for which experimental cross sections are known for a wide range of 
the incident energy, from the thermal neutron energy up to the pion-production threshold. 
We therefore calculate here both u — d reaction and np -^ 'yd cross sections in the same 
formalism and use the latter to gauge (at least partially) the reliability of our model. 

III. CALCULATIONAL METHODS 

A. Multipole expansion of hadron current 

To evaluate the two-nucleon matrix element of the hadron current, we first separate the 
center-of-mass and relative wave functions, 

< ri,r2 I d{P) >= e^^-^^d(r), < r,,r2 \ NN{P') >= e^^'-^^p,(r), (35) 

where r = ri — r2 and R = ^^^~^, and ipd and ipp' represent, respectively, the deuteron wave 
function and a scattering-state wave function with asymptotic relative momentum p'. Then 
the matrix element of the hadron current for charged-current reaction is given by 



jf^ =< NN{Py I jf^(O) I d{P) > = fdr ^;,(r) f dR e-^'^-^jf^(O) 



Mr). (36) 



As for the neutral-current reaction, we just replace Jf^*" with Jj^*" . In the following equations, 
Jx without superscript applies for both NC and CC. Eliminating the dependence of the 
current Jx{x) on the center-of-mass coordinate, R, we can write 

Jx =< MJd^ e^''-"JA(a^)|^d >, (37) 

where JA(aj) = J\{x)\ii^o. Similarly, we define Va(x) = VA(a;)|H=o, and^A(a;) = Ax{x)\ii^Q. 
We now introduce the standard multipole expansion of the nuclear currents WM . The mul- 
tipole operator for the time component of a current is defined by 

13 



Tt\J) = J dxjj{qx)YjM{x)Mx), (38) 

where jj{qx) is the spherical Bessel function of order J, q = \q\, and x = x/\x\. The electric 
and magnetic multipole operators are defined by 

Ti^'iJ) = - fdxVx [jj{qx)YjjM{x)] ■ J{x), (39) 

q J 

Til\j) = Jdxjj{qx)YjjM{x) ■ J{x), (40) 

where Yjlm{x) are vector spherical harmonics. The longitudinal multipole operator is 
defined by 

Tl''{J) = - f dxV[jj{qx)YjM{x)] ■ J{x). (41) 

q J 

Using the conservation of the vector current, the longitudinal multipole operator of the 
vector current can be related to the charge density operator as 

T/°(V) = --T^°(V). (42) 

An explicit form of the electric multipole operator for the vector current is given by 



Ti^iy) = -i^,^^ j dxjj+,{qx)Yjj+,M{&) ■ V{x) 

+'V ^J^ J dxjj-i{qx)Yjj^iM{x) ■ V{x). (43) 

Here again we can use the current conservation to rewrite Eg .(^31) into a form that has the 
correct long wavelength limit of an electric multipole operator: 



TnV) = -^|■^-^n'\V)-^f-^ j dx3j^,{qx)Yjj^,M{&)-V{x). (44) 



B. Cross sections 

As explained earlier, we calculate the cross sections for v /v{k) + d{P) -^ l{k') + Ni{p[) + 
N2{p'2) in the laboratory system. Following the standard procedure, we obtain the cross 
section for the CC reaction as 

14 



^^ = E ^ 7^, 2 ^ ' '^ \l 3x^?dk'dp\dp'2, (45) 

and the cross section for the NC reaction as 

da = Y. ^-^^^^^f^-^^^rjn'dk'dp\dp',. (46) 

The matrix elements, l^ and jx, have been defined in Eq.(0) and in Eg . (|36|) , respectively. 
In Eg. (|45D , we have included the Fermi function F{Z,E'g) [|^ to take into account the 
Coulomb interaction between the electron and the nucleons. In fact, this factor is relevant 
only to the Ue + d^e^+p + p reaction, for which we should use F{Z = 2, E'^); for the 
Ue + d~^e^ + n + n reaction we have F{Z = 0, E'^) = 1. 

Substitution of the multipole operators defined in Eg.(^) ~Eg.(^l]) leads to 



I'jx = E 47rz^°(-l)^^° 

Jo Mo 

X < V^p/ I ^T^oMo^J^o-AIo ^ rpJoMo^Jo~Mo ^ J.JoMo^Jo-Mo ^ rpJoMo^Jo-Mo^ | ^^ ^ ^ ^4^) 

where the lepton matrix elements are given as 

ii"" = YjMiq) f, (48) 



?if = YjjMiq) ■ Z, (50) 



To proceed, we use a scattering wave function of the following form: 

i^p'{r)= Yl 47r(l/2,si,l/2,S2|5/i) (l/2,ri,l/2,r2|T,T,) (Lm5/i|JM) 

L,S,J,T 

X t'-Yl^ip') ^PlsjAt) (52) 

with 

i^LsMr) = ^ ~ ^~^l ^ ^ J2 yL'sAr) Ri,L;s{r) Vt,t., (53) 

yLSj{r) = [YL{r)®XsU, (54) 
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where xs (jIt) is the two- nucleoli spin (isospin) wave function with total spin S (isospin T). 
The above wave function is normalized in such a manner that, in the plane wave limit, it 
satisfies 

Ri',L,si'')^jL{p'r)5L,v. (55) 

The partial wave expansion of the scattering wave function [Eq.(^2|)] gives 

^'JA = E E (-1)"'° ^'°"^^ft=T( V2, ^1, 1/2, s,\S^,) (1/2, n, 1/2, r^lT, T.) 

x(lmdJoMo|JM)(Lm5/x|JM) Fi,^(p') E < T^" > 4°"^^°, (56) 

X=C,£,L,M 

where tti^ is the z-component of the deuteron angular momentum. We have used here a 
simplified notation 

<0^°> = <^LSJT\\0'''\\iJd> (57) 

for the reduced matrix element defined by 

<J'M'\ 0-^°*^° \JM> = , ^ (J,M,Jo,Mo\J',M') < J'||0^°|| J>, (58) 

where 0'^°'^° are the multipole operators that appear in Eq.(^)~(^T]). 



1. Cross sections for charged- current reaction 

For the CC reaction, observables of interest are the total cross section and the lepton 
differential cross sections. We therefore integrate Eq. ( ^5]) over the momenta of the final two 
nucleons. The evaluation of the phase space integrals and the relevant kinematics are briefly 
described in the Appendix. According to the Appendix, Eq.(^) leads to 

da = ^"""l F{Z, E',) |M|25(Md + k-E[- P">)J p'^dp'k'^dk'dQk', (59) 

where 
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|M|2 = J2 { I < ^C°(V) > |'(1 + k-l3+'^il-k-(3 + 2q-(3q-k)- — q-{k + f3)) 

LSJ,Jo 1 1 

+ 1 < T^-iA) > 1^(1 + fc ■ /3) + I < T/°(^) > |2(1 -k-f3 + 2q-f3q-k) 

+2Re[< T^^iA) >< Tl°{A) >*] q ■ {k + f3) 

+ [\ < Ti!{V) > P + I < Ti°(V) > P + I < T//(^) > P + I < Ti°(^) > |2] (1 - q . fc q . /3) 

T2Re[< Ti^{V) >< Ti°{A) >* + < T//(^) x Ti°(V) >*] g • (fc - /3) } . (60) 

In the above, k' = \k'\ and /3 = k'/E'f, p' is the relative momentum of the final two 
nucleons, and p' = \p'\. Of the double sign in the last line of Eg. (|60|) , the upper (lower) 
sign corresponds to the z/ {u) reaction. The appearance of the factor J in Eq.(^) needs an 
explanation. As discussed in the Appendix, when relativistic kinematics is adopted, there 
arises a Jacobian, J, associated with the introduction of p' but it is a good approximation 
to use J, the angle-averaged value of J. 

For the total cross section, the use of relativistic kinematics gives 

J J Stt l + E'^{l-kcoseL/k')/JPI^^ + q^ 

where T is the kinetic energy of the final NN relative motion and 6l is the lepton scattering 
angle (cos 6^ = k- k') in the laboratory frame. If instead we use non-relativistic kinematics, 
the results would be 

J J ^ ' Stt l + E'^{l-kcoseL/k')/{Mm + MN2) ^ ' ^^i ' ' ^ ^ 

where M^i is the mass of the i-th nucleon, and M^ is the reduced mass of the final NN 
system. 

Eq. (|59|) also leads to the double differential cross sections for the z/g + ci— >-e^+]9 + p 
reaction: 

_|^ ^ ^-^^FiZ, E',)m^E',41fTV^ W?. (63) 

The electron energy spectrum and the electron angular distribution are obtained from 
Eq.(|63l) as 
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da r ,^ ( d'^cr \ da f ,„, f d 



a 



d^'^' h?^-T^ :J?^= dE'J-—- _. (64) 



dEg J \dQk'dEfJ^,^^. dQy J \d^k'dE[j ^ 



2. Cross sections for neutral- current reaction 

The total cross section for the NC reaction can be calculated in essentially the same 

manner as above. The result is 

r r 2Gl JE'(.fPJ^/2]p'k' 

a= dT dicose^)^^ iS^J^LJJL^^=\M\\ (65) 

J J ' ' Svr 1 + E',{1 - k cos dL/k>)/^P';^ + q' 

where |Mp is given by Eq.(|60D with, however, the charged current replaced by the neutral 
current. By contrast, in calculating neutron differential cross sections we can no longer 
integrate over the relative momentum of the final nucleons. We therefore work with the 
following expressions: 

;/^=/^^-#^#^^-T% E l^'^^'l'' (66) 

dUj^JTn J 3(27r)5 Ep - p; ■ fc' ^,X«p 

where we have indicated explicitly averaging over the initial spin and summing over the final 
spins. The energy and momentum of the final proton (neutron) are denoted by {E'^,p'a) 
with a = p {a = n); Tn is the kinetic energy of the neutron. The neutron energy spectrum 
and the neutron angular distribution are then evaluated as 

^^ ' dn,^ ( ^'" ] -^= fdT^ ( ^'" ] . (67) 



din J V'^^^Pn'^-'"/ Eq.(66) '^^ ^Pn •' V*^^ ^Pn -^ " / Eq.@ 

The calculation of the total cross section for the np -^ 'yd reaction follows essentially the 
same pattern as that of the u — d total cross section, and therefore we forgo its description. 

IV. NUMERICAL RESULTS 

A. Radiative capture of neutron on proton 

To test the nuclear currents and wave functions used, we first discuss the capture rate 
for np —>■ 'yd. Thermal neutron capture is a well known case for testing exchange currents 



T9| , |20[] . This reaction is dominated by the isovector magnetic dipole transition from the 



^5*0 np scattering state. With the use of the ANLV18 potential, our PhLA calculation gives 
a{np — > '-fd) = 335.1 mb, with both the lA and MEX currents included. This is in good 



agreement with the experimental value, a{np — > ^dY^'^ = 334.2 ± 0.5 mb ||^. With the lA 
contribution alone, our result would be a{np -^ 'yd)^'^ = 304.5 mb. The 10% contribution 
of the exchange current is due to the pion, pair and A-currents. 

Going beyond the thermal neutron energy regime, we give in Fig.|l|the calculated a{np — >■ 
'yd) as a function of the incident neutron kinetic energy, T„. The experimental data in Fig.|l| 
have been obtained from either the neutron capture reaction itself |^^ or its inverse process 



48| , |49[| , using detailed balance for the latter. We can see that our results describe very 



well the energy dependence of a{np -^ 'ydy^^ all the way up to T„ k, 100 MeV. The figure 
indicates that the electric dipole amplitude starts to become important around T„=100 
keV. In the higher energy region we should expect deviations from the long-wave length 
limit of the electric dipole operator, and therefore the good agreement of our results with 
the data suggests that the description of the electric multipole is also satisfactory.]] The fact 
that our PhLA calculation with no ad hoc adjustment of the input parameters is capable of 
reproducing a{np -^ 'ydy^^ for a very wide range of the incident energy gives us a reasonable 
degree of confidence in the basic idea of PhLA and the input parameters used.0 Of course, 
strictly speaking, the electromagnetic and weak-interaction processes do not probe exactly 
the same sectors of PhLA, but the remarkable success with a{np -^ 'yd) gives, at least, 
partial justification of our PhLA as applied to weak-interaction reactions. Noting that the 
dominant axial MEX current due to A-excitation is related to the A-excitation MEX current 
for the vector current (we need only replace (/y + /A/)/2MAr with /^ ), we evaluate the former 



^Since our treatment here does not include pion production, our results should be taken with 
caution above the pion production threshold. 



Another similar success of PhLA is known in the d{e, e')np reaction |2C|. 
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with the same input parameters as used in calculating a{np -^ •yd). 

B. Cross sections o{ v — d reactions 

We now present our numerical results for the z/(z/) — d reactions. In what follows, the 
"standard run" represents our full calculation with the following features. The ANLV18 
potential |^ is used to generate the initial and final two-nucleon states and the final two- 
nucleon partial waves are included up to J = 6. For the transition operators, we use the 
lA and MEX operators described in Sec. |I|; the Siegert theorem is invoked for the electric 
part of the vector current. As regards the single- nucleon weak- interaction form factors, we 
employ the most updated parameterization given in Eqs.(^)-(^). The final two-nucleon 
system is treated relativistically in the sense explained in Appendix.^ Our numerical results 
will be given primarily for our standard run; other cases are presented mostly in the context 
of examining the model dependence. 

1. Total cross sections for v — d and D — d reactions 

We give in TABLE H and Fig.0 the total cross sections, obtained in our standard run, for 
the four reactions: Ued — > e~pp, u^d -^ t^x^p, Ugd — > e~^nn, u^d -^ Vx^p. The cross sections 
are given as functions of E^, the incident v/v energy, from the threshold to E^ = 170 MeV. 
It should be mentioned that towards the highest end of Ei, considered here, pion production 
sets in but the present calculation does not include it. 

It is informative to decompose the total cross section into partial-wave contributions. 



^ We must emphasize that our calculation takes account of "relativity" only in certain aspects of 
kinematics. Going beyond this is out of the scope of this paper. 
^^The numerical results reported in this article are available in tabular and graphical forms at the 



website: <[http://nuc003.psc.sc.edu/rkubodera/NU-D-NSGK>. 
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TABLE H shows the relative importance of the two lowest partial waves in the final two- 
nucleon state; denoting the contributions to the total cross section from the ^5*0 and ^Pj 
states by ai^So) and I]jO"('^-Pj)) respectively, we give in TABLE ^ the ratios, cr("'^S'o)/cr(all) 
and Z]j=o '^('^-Pj)/c"(8'll)) as functions oi E^. Here (j(all) denotes the sum of the contributions 
of all the partial waves; in fact, it is sufficient to include up to J = 6 even for E^ = 170 
MeV, where the summed contribution of higher partial waves (J > 6) is found to be less than 
1%. The table reconfirms that, in the low-energy region, the Gamow- Teller (GT) amplitude 
due to the ^5*0 final state gives a dominant contribution. It is therefore important to take 
into account the A-excitation axial- vector current, which gives a main correction to the lA 
current. As mentioned, in our approach, the coupling constant determining the A-excitation 
MEX current is controlled by the np -^ 'yd amplitude. As E^ increases, the ^Pj final states 
become as important as the ^5*0 state, and therefore 1^ type multipole operators arising 
from the vector as well as axial-vector currents start to play a significant role. In this sense 
it is reassuring that the validity of our model for the electric dipole matrix element in this 
energy region has been tested in the photo-reaction. 

Turning now to TABLE ^, we give in the second column labeled "lA" , the ratio of the 
total cross section obtained with the use of the lA terms alone to that of our standard run. 
We see that, at the low energies, the MEX contribution is about 5% of the lA contribution. 
As Ej, increases, the relative importance of the MEX current contribution augments and it 
can reach as much as 8% in the high energy region. The third column {+Amex) in Table 



rT| gives the cross section that includes the contribution of the space component of the axial 
exchange current, while the fourth column {+AKDRfl) gives the results that contain the 
additional contribution of the time component of the axial exchange current. It is clear that 
the MEX effects are dominated by +Amex] the axial-charge contribution is very small for 
the entire energy range considered here. The last column {+V'mex) i'^ Table |T| gives results 
obtained with the use of the full vector exchange currents, Eq.(^), i.e., without invoking 
the Siegert theorem. The numerical difference between the two cases (with or without the 
Siegert theorem imposed) is found to be very small; the difference is practically zero for lower 
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values of Ei, and, even at the higher end of E,^, it is less than 1%. Thus the Siegert theorem 
allows us to take into account implicitly most part of the MEX for the vector current .[^ 

In order to compare our cross sections with those of the previous work, we give in 
TABLE |rV| the ratios of the cross sections reported in YHH [|T^ and in KN [|12[ to those 



of our standard run; the second column gives a{Y]:M.) / a{standard run), while the third 
column shows a{KN) / a{standard run). In the solar neutrino energy region, one can see 
that the results of our standard run agree with those of KN |jl2| within 1% except for the 



Ued -^ e^pp reaction near threshold, wherein the discrepancy can reach 2%. As the incident 
energy becomes higher, our results start to be somewhat larger than those of KN, and the 
difference becomes about 6% towards the higher end of E^. This variance arises largely 
from the cutoff mass in the form factor G'a(q'^), which accounts for 3-4% difference.0 The 
remaining ~2% difference is due to our use of relativistic kinematics and the inclusion of the 
contributions from higher partial waves and from the isoscalar current which were ignored 
in the previous study. We have done an additional calculation by running our code adopting 
the same approximations and the same input parameters as in KN, and confirmed that the 
results agree with those of KN within 1% in the high energy region as well.[^ 



^^ In our approach, which uses phenomenological nuclear potentials, the conservation of the vector 
current is not strictly satisfied. A measure of the effect of current non-conservation may be provided 
by comparing two calculations, one with the Siegert theorem implemented and the other without. 



The results in Table III indicate that numerical consequences of the current non-conservation are 
practically negligible in our case. 

^^The value of the cutoff mass, niA, in [^ was deduced from an experiment involving a deuteron 
target and therefore it may involve nuclear effects. It seems worthwhile to reanalyze the data 
taking into account possible nuclear effects. Another potentially useful source of information on 



niA is low-energy pion electroproduction |50|. 



^^The precision of our numerical computation of the cross sections is also 1%. 
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On the other hand, the cross sections of YHH |]T0| are about 5% smaller than those of 
our standard run even at the low energy. This reflects the fact that YHH did not include 
the MEX contributions (except for the term that could be incorporated via the extended 
Siegert theorem). Indeed, comparison of the YHH cross sections with the entries in the 
second column labeled "lA" in TABLE |ITT| indicates that, if we drop the explicit MEX 
terms in our calculation, the resulting cross sections in the solar energy region agree with 
those of YHH within ~1%. 

We next consider the NN-potential dependence of the cross sections. The fourth column 
labeled "Reid93" in TABLE |I^ gives the ratio of the total cross section obtained with the 
use of the Reid93 potential ||4^ to that of our standard run; the fifth column gives a similar 
ratio for the case of the NIJ II potential [^. We note that the dependence on the nuclear 
potentials is within 1% for all the reactions and for the entire energy region under study. 
Since all the potentials used here describe the NN scattering data to a satisfactory degree, 
it is probably not extremely surprising that all these modern realistic NN potentials give 
essentially identical results for v — d cross sections, but the present explicit confirmation is 
reassuring. 

In our calculation the strength of the A-excitation exchange current, which contributes 
both to the Gamow- Teller and Ml transitions, is monitored by the empirical values for 



a{np — > 7(i). Meanwhile, Carlson et al. |2J], in estimating the solar pp-fusion cross section, 
used the tritium /?-decay rate to fine-tune the vrA^A coupling constant that features in the 
Gamow- Teller exchange current. This method turns out to yield somewhat "quenched" A- 
excitation MEX effects in the pp fusion. It is therefore of interest to study the consequences 
of this second method for the v — d reactions. In the last column labeled "A(CRSW)" of 



^^There is 2% variance for the Ved -^ e^nn cross section near threshold (not shown here); this 
is however very likely to be attributable to the fact that the n-n scattering length is not exactly 
reproduced by the potentials other than ANLV18. 
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TABLE 1^, we give the ratio of the cross sections obtained with the use of the A-current 



employed in [Q to those of our standard run. In the solar energy region this ratio is found 
to be 0.96 - 0.97, or the MEX contribution relative to the lA term is 2%, instead of 5% 
found in our standard run. This reduction is primarily due to the smaller ttNA coupling 



constant in [^. At higher neutrino energies, the use of the the A-current employed in ||2^ 
leads to a ~4% MEX effect relative to the lA term, to be compared with the ~8% effect 
found in our standard run. Thus, in general, if we adopt the approach taken in [^, the 



importance of the MEX effect relative to the lA contribution will be reduced by a factor of 
~2 as compared with the result of our standard run. 



As emphasized by Bahcall et al. ||3^ , one of the crucial quantities in neutrino oscillation 
studies at the SNO is the double ratio [NC]/[CC], where [NC] ([CC]) itself is the ratio 
of the observed neutrino absorption rate to the standard theoretical estimate for the NC 
(CC) reaction rate. This implies that the reliability of theoretical estimates for the ratio 
R = a{NC)/a{CC) = aiyd —>■ vnp) / a{Vf.d -^ e~pp) is extremely important. We give in 
Table |V| the values of R resulting from the various models considered in this paper. Since 
our primary interest here is to examine the model dependence of R, we choose, in Table |V], 
to normalize R by R standard rum the value corresponding to our standard run] Rstandard run 
itself is shown in the second column of the table. We learn from Table that all the models 
studied give essentially the same R; deviations from Rstandard run are at most ~1%. Thus, 
the largest source of model dependence in our work due to the A-exchange current cancels 
out by taking the ratio between the NC and CC reactions. 

2. Differential cross sections for ttie electron 

We now discuss three types of electron differential cross sections for the Ue+d —>■ e^-\-p+p 
reaction: (i) the energy spectrum, da/dE'^ in Eg. (|6^) , (ii) the electron angular distribution, 
da /dVty in Eg. (p^) , and (iii) the electron double differential cross sections, d'^a /dE'^dVty in 
Eg. (^3]). Although this kind of information must be implicitly contained in the computer 
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codes used in the existing work |]8|- |T^] , its explicit tabulation has been lacking in the liter- 
ature. It seems very useful to make these differential cross sections readily available to our 
research community. However, a trivial but nonetheless serious problem is that the required 
amount of tabulation is enormous. We therefore present here some representative results, 
relegating the bulk of tabulation to a web site.[^ For four values of the incident neutrino 
energies, E^ = 5, 10, 20 and 150 MeV, we give the electron-energy spectra, da/dE'^, in 
Fig.|| and the electron angular distribution, da/dfik', in Fig.§. We note that the electron 
spectrum in Fig.|^ exhibits a "cusp-like" structure for £',^=150 MeV. This feature, which is in 
fact common for Ei, > 100 MeV, probably calls for an explanation. For a given value of E^y, 
we can separate the electron energy E'^ into two ranges: E'^ < E'f or E'^ > E'f, where E'f is 
the point above which the electron scattering angle 6l cannot any longer cover the full range 
[0, n] for a kinematic reason.]^ The "cusp-like" structure occurs at E'^ = E'f due to interplay 
between the change in the range in the phase space integral and the momentum dependence 
in the transition matrix element for the final ^5*0 channel. This structure, however, is not 
a cusp in the mathematical sense. Enlarging the scale of the abscissa, we can confirm that 
the actual curve is a rapidly changing but non-singular one. It turns out that for higher 
values of E^ we need more scale enlargement before the curve starts looking smooth to the 
eye. This is the reason why, for a fixed abscissa scale (as adopted in our illustration), the 
case corresponding to the high incident energy tends to exhibit more "cusp-like" behavior. 
Regarding the electron angular distribution (Fig.^, we note that at low neutrino energies 
the electrons are emitted in the backward direction, carrying most of the available energy. 
The angular distributions for the lower incident energies are reminiscent of that for a Gamow- 
Teller /?-decay between two bound states. If we simplify the expression for the electron 
differential cross section, (Eq.(|63D), by dropping all the partial waves other than ^5*0 and by 

^5 See footnote |l^. 
^^See the Appendix. 
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retaining only the leading-order Ganiow- Teller matrix element, then we have 

/^2 2 n 

da ~ 1T3 "^ flMpp'k'^F{Z, E'^) (3 - P cos 9 L)I^dk'dnk', (68) 

where / is the relevant radial integral. Since /? ~ 1 and F ~ 1, if we tentatively treat / as 
a constant, we have a simple expression 

da oc p'k'^ (3 - cos eL)dk'dVLk' . (69) 

In fact, the electron angular distributions for low incident neutrino energies can be simulated 
to high accuracy by Eg. (p9|) ; see the dotted lines in Fig.||. Thus, although the radial integral 
I may in fact depend strongly on the kinetic energy of the NN relative motion, the numerical 
results for da /dVLy at the low energies can be conveniently simulated by the simple phase- 
space formula, Eq.(|69D. 

As for the electron double differential cross sections, d'^a/dE'^dQ^', Eq.(^), even pre- 
senting some typical cases is impractical because of the bulkiness of the tables. We therefore 
relegate their tabulation completely to the web site the address of which is given in footnote 



10. 



3. Neutron energy spectrum and angular distribution 

Finally, we consider the neutron energy spectrum, da/dTn and the neutron angular 
distribution, da/dQn in Eg. (|67D , for the u + d ^ u + p + n reaction. For E^ = 5, 10, 20, and 
50 MeV, we show da/dT^ in Fig.|^ and da/dfln in Fig.|[ Once again, we relegate a complete 
tabulation of our numerical results to the web site mentioned in footnote 0. We see from 
Figs.^ and ^ that the neutron energy spectrum has a peak near the lower end, and that, 
unlike the electrons, the neutrons are emitted in the forward direction. 

V. DISCUSSION AND SUMMARY 

Based on a phenomenological Lagrangian approach (PhLA), we have carried out a de- 
tailed study of the u — d reactions and provided the total cross sections and the differential 
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cross sections for the electrons and neutrons, from threshold to E,y = 170 MeV. We have 
examined the influence of changes in various inputs that feature in our PhLA. In particular, 
we have studied to what extent the use of the modern NN potentials affects the results. We 
have also examined the influence of the use of the updated input concerning the nucleon 
weak-interaction form factors. The vertex strength that governs the A-excitation axial- 
vector exchange current has been monitored using the photo-reaction. We have also studied 
the consequence of employing the vertex strength determined with the use of the tritium 
/3-decay strength p^ . 

For the solar energy region, E,y < 20 MeV, the results are summarized as follows. By 
comparing our new results with those in the literature, we have confirmed that the total 
V — d cross sections are stable within 1% precision against any changes in the input that 
have been studied, except for somewhat higher sensitivity to the strength of the A-current 
(see below). The same stability should also exist for the differential cross sections described 
in this paper. The MEX axial-vector current in our standard run increases the total cross 
sections by ~5% from the lA values; we have used the np — > ^d reaction to monitor the 
dominant part of our MEX current. Meanwhile, Carlson et al. |Q, in estimating the solar 



pp-fusion cross section, used the tritium /3-decay lifetime to monitor a vertex strength that 



features in the Gamow- Teller exchange current. The results of pj] indicates that adjusting 
the MEX strength using the tritium /5-decay rate could lead to a somewhat reduced MEX 
amplitude. If we use the A-excitation axial current renormalized by the tritium /3-decay 
m, the MEX current correction to the lA term, [cr (lA + MEX) - (T(IA)]/a-(IA), turns out 
to be ~2%, instead of 5% as in our standard run; see the column labeled A(CRSW) in 
TABLE 1^. The difference between our standard run and A(CRSW) represents the range 
of uncertainty in the present PhLA calculation. We therefore consider it reasonable to use, 
as the best estimates of the low-energy vd cross sections, the values given by our standard 
run and attach to them a possible overall reduction factor k, with k ranging from 0.96 to 
1. In this language, the "Icr" uncertainty adopted by Bahcall et al. |^3[ corresponds to 
K = 0.95 — 1, which represents the difference between the cross sections given in YHH [l^ 
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and KN [0. We have shown that in the ratio R = o"(NC)/o"(CC) the model dependence is 
reduced down to the 1% level (see TABLE [V|). 

At higher incident neutrino energies, the results obtained in our standard run are some- 
what larger than those of KN, and the difference reaches ~6% towards E,y=150 MeV. This 
difference is caused largely by the updated value for the axial-vector mass. The effect of 
relativistic kinematics, as discussed here, has ~1% effect on the cross sections. The contri- 
butions of the isoscalar current, which so far has been totally ignored in the literature, is 
found to be of 1% even at E'j^ ~ 150 MeV. The importance of the MEX currents relative to 
the lA contributions increases monotonically as E^, augments. Towards £',^=150 MeV, the 
MEX to lA ratio, [cr(IA + MEX) - a(IA)]/o-(IA), reaches ~8% in our standard run while 
this ratio is ~4% in the case of A(CRSW). 

As mentioned earlier in the text, the numerical results of this work are fully documented 



in tabular or graphical form at the website referred to in footnote [T0|. It is hoped that those 
tables and graphs are of value for the ongoing and future neutrino experiments that involve 
deuteron targets. 
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APPENDIX A: PHASE SPACE INTEGRAL AND KINEMATICS 



We briefly explain the derivation of the cross section formula Eq.(p9D starting from 



Eq.(|5). The phase space integral in Eg. (p9|) is 

/ = S\k + P-k' -P') dp\dp'2dk' 



= 5{E, + Md-E[- ^P'2 + p;2 ) dp'^dk', (Al) 

where p' l = (p'l — p'2)/2 and P' = q = k — k' . 

The scattering energy of the final NN distorted wave is given by their center-of-mass 
energy Wnn = \ Pn- The relative momentum in the center-of-mass system, p"^, is given by 



Lorenz-transforming the relative momentum in the lab system as [^| 

p" = K P'l- (A2) 

The magnitude of p' is related to Wnn as 



Wnn = Vp" + Ml^ + Vp'' + M2,2, (A3) 

where M^i is the mass of the z-th nucleon in the final state. The integral over the momentum 



p' L is then replaced by integration over p', which gives rise to a Jacobian [^ , 

dp'L = Jdp', (A4) 

with 



Wmn{E[ + E',) 



J = ..r ,^, ' ^.^ ^ (A5) 



where E[ is the energy of the i-th nucleon in the lab system. Although J depends on the 
direction of p', we approximate it by J = ^ / JdQp> ; through a plane- wave calculation, we 
have confirmed that this is a good approximation in the energy region of our concern. The 
phase space integral is then given as 



/ = 6{E, + Md-E[- ^P'2 + p;2 ) jdp'dk', (A6) 
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which leads to Eq.(p9D. 

The kinematically allowed domain of the integral dp'dk' is determined by a standard 

procedure. We give here the results for the electron energy spectrum, Eq.(0), for the 

Ve + d^>-e~+p + p reaction. The threshold neutrino energy E*J^ for this reaction is given 

by 

, (2Mp + M, + me)(2Mp - M, + m,) 
^^ = 2M, ■ ^^^^ 

We may specify the allowed region of the electron energy E'^ by giving the conditions on the 

electron momentum k'\ these conditions are 



(A8) 



where 



< A;' < A;^ for E^ > E^ 

k'_ <k'<k'^ for El>E^> El'', 

{2Mp + Ma-m,){2Mp-Md + m,) 
""' = 2(M, - m.) ' (^') 



and 



k'. = ^ . (AlO) 

with W'^ = {P + k)l and X = Mj + 2E^Md - AM^ + m^. For given values of E^ and E'^, 

the electron scattering angle 6^ is restricted as 

f 2E'(Md + E^)~X] ^ ,, , 

max|-l, ^ } - '°'^^ - ^' ^^^^^ 

and the NN scattering energy is specified by p' given as 



p' = -^X + 2E,k' cos Ol - 2E'^{Md + E,). (A12) 

For Ej^ = 150 MeV, the allowed ranges of cos di and p' are plotted as functions of E'^ in 
Fig. 7 and Fig. 8, respectively; the dotted area in each figure represents the allowed region. 
At -Eg = E'^, the constraint on 6l sets in and the minimum value of p' becomes zero. E'^ is 
determined from the condition: 

2E'f{Md + E,) + 2E,k'' -X = 0. (A13) 
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TABLES 
TABLE L Total cross sections for u — d reactions in units of cm^. The "— x" in parentheses 
denotes 10"''; thus an entry hke 4.279 (-47) stands for 4.279 xlO^^'^ cm^. 



E^ [MeV] 


v d ^ 


J p n 


ly d ^ 


y p n 


Vf, d ^ t 


3 pp 


V(, d ^ e'^ n n 


2.0 


0.000 


(0) 


0.000 


(0) 


3.603 


;-45) 


0.000 ( 0) 


2.2 


0.000 


(0) 


0.000 


(0) 


7.833 


;-45) 


0.000 ( 0) 


2.4 


4.279 


:-47) 


4.248 


-47) 


1.404 


;-44) 


0.000 ( 0) 


2.6 


4.258 


;-46) 


4.222 


-46) 


2.242 


;-44) 


0.000 ( 0) 


2.8 


L457 


;-45) 


1.443 


-45) 


3.315 


;-44) 


0.000 ( 0) 


3.0 


3.355 


;-45) 


3.320 


-45) 


4.639 


;-44) 


0.000 ( 0) 


3.2 


6.286 


;-45) 


6.213 


-45) 


6.228 


;-44) 


0.000 ( 0) 


3.4 


L038 


;-44) 


1.025 


-44) 


8.095 


;-44) 


0.000 ( 0) 


3.6 


L574 


;-44) 


1.553 


-44) 


1.025 


:-43) 


0.000 ( 0) 


3.8 


2.246 


;-44) 


2.213 


-44) 


1.271 


;-43) 


0.000 ( 0) 


4.0 


3.060 


;-44) 


3.012 


-44) 


1.547 


:-43) 


0.000 ( 0) 


4.2 


4.024 


;-44) 


3.956 


-44) 


1.855 


;-43) 


1.115 (-45) 


4.4 


5.142 


;-44) 


5.049 


-44) 


2.196 


:-43) 


4.554 (-45) 


4.6 


6.420 


;-44) 


6.297 


-44) 


2.570 


:-43) 


1.010 (-44) 


4.8 


7.860 


;-44) 


7.702 


-44) 


2.978 


;-43) 


1.787 (-44) 


5.0 


9.468 


;-44) 


9.267 


-44) 


3.420 


;-43) 


2.799 (-44) 


5.2 


1.125 


;-43) 


1.100 


-43) 


3.897 


;-43) 


4.059 (-44) 


5.4 


1.320 


:-43) 


1.289 


-43) 


4.410 


:-43) 


5.578 (-44) 


5.6 


1.533 


;-43) 


1.495 


-43) 


4.959 


:-43) 


7.364 (-44) 


5.8 


1.763 


;-43) 


1.718 


-43) 


5.544 


;-43) 


9.427 (-44) 


6.0 


2.012 


;-43) 


1.958 


-43) 


6.166 


;-43) 


1.177 (-43) 


6.2 


2.279 


;-43) 


2.215 


-43) 


6.825 


;-43) 


1.441 (-43) 


6.4 


2.564 


:-43) 


2.490 


-43) 


7.522 


;-43) 


1.733 (-43) 


6.6 


2.868 


;-43) 


2.782 


-43) 


8.258 


;-43) 


2.056 (-43) 


6.8 


3.191 


;-43) 


3.092 


-43) 


9.031 


;-43) 


2.409 (-43) 


7.0 


3.532 


;-43) 


3.419 


-43) 


9.843 


;-43) 


2.792 (-43) 


7.2 


3.893 


:-43) 


3.764 


-43) 


1.069 


;-42) 


3.206 (-43) 


7.4 


4.273 


:-43) 


4.126 


-43) 


1.159 


:-42) 


3.652 (-43) 


7.6 


4.672 


;-43) 


4.506 


-43) 


1.252 


;-42) 


4.127 (-43) 


7.8 


5.091 


;-43) 


4.904 


-43) 


1.349 


;-42) 


4.635 (-43) 


8.0 


5.529 


:-43) 


5.320 


-43) 


1.450 


;-42) 


5.175 (-43) 


8.2 


5.987 


;-43) 


5.754 


-43) 


1.555 


;-42) 


5.746 (-43) 


8.4 


6.464 


;-43) 


6.206 


-43) 


1.664 


:-42) 


6.349 (-43) 


8.6 


6.961 


;-43) 


6.676 


-43) 


1.777 


;-42) 


6.984 (-43) 


8.8 


7.479 


;-43) 


7.163 


-43) 


1.894 


;-42) 


7.652 (-43) 


9.0 


8.016 


:-43) 


7.669 


-43) 


2.016 


;-42) 


8.351 (-43) 


9.2 


8.573 


;-43) 


8.193 


-43) 


2.141 


;-42) 


9.082 (-43) 


9.4 


9.150 


;-43) 


8.735 


-43) 


2.271 


:-42) 


9.846 (-43) 
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TABLE I. (continued) Total cross sections for v — d reactions in units of cm^. The 



-X m 



parentheses denotes 10 ^; thus an entry hke 4.279 (-47) stands for 4.279 xlO 



-47 pj^2_ 



Eu [MeV] 



V d ^ V p n 



u d ^ u p n 



Ve d ^ e p p 



Ve. d 



n n 



9.6 

9.8 

10.0 

10.2 

10.4 

10.6 

10.8 

11.0 

11.2 

11.4 

11.6 

11.8 

12.0 

12.2 

12.4 

12.6 

12.8 

13.0 

13.5 

14.0 

14.5 

15.0 

15.5 

16.0 

16.5 

17.0 

17.5 

18.0 

18.5 

19.0 

19.5 

20.0 

20.5 

21.0 

21.5 

22.0 

22.5 

23.0 

23.5 

24.0 

24.5 



9.747 


'-43) 


1.036 


'-42) 


1.100 


'-42) 


1.166 


'-42) 


1.234 


'-42) 


1.304 


'-42) 


1.376 


'-42) 


1.450 


'-42) 


1.526 


'-42) 


1.604 


'-42) 


1.684 


'-42) 


1.767 


'-42) 


1.851 


'-42) 


1.938 


'-42) 


2.026 


'-42) 


2.117 


'-42) 


2.210 


'-42) 


2.305 


'-42) 


2.551 


'-42) 


2.811 


'-42) 


3.084 


'-42) 


3.371 


'-42) 


3.671 


'-42) 


3.984 


'-42) 


4.311 


'-42) 


4.651 


'-42) 


5.006 


'-42) 


5.374 


'-42) 


5.755 


'-42) 


6.151 


'-42) 


6.560 


'-42) 


6.984 


'-42) 


7.421 


'-42) 


7.872 


'-42) 


8.338 


'-42) 


8.817 


'-42) 


9.311 


'-42) 


9.819 


'-42) 


1.034 


'-41) 


1.088 


'-41) 


1.143 


:-41) 



9.294 


-43) 


9.872 


-43) 


1.047 


-42) 


1.108 


-42) 


1.171 


-42) 


1.236 


-42) 


1.303 


-42) 


1.372 


-42) 


1.442 


-42) 


1.514 


-42) 


1.588 


-42) 


1.664 


-42) 


1.741 


-42) 


1.821 


-42) 


1.902 


-42) 


1.985 


-42) 


2.069 


-42) 


2.156 


-42) 


2.379 


-42) 


2.614 


-42) 


2.860 


-42) 


3.117 


-42) 


3.385 


-42) 


3.663 


-42) 


3.953 


-42) 


4.253 


-42) 


4.564 


-42) 


4.886 


-42) 


5.218 


-42) 


5.561 


-42) 


5.915 


-42) 


6.279 


-42) 


6.653 


-42) 


7.038 


-42) 


7.434 


-42) 


7.839 


-42) 


8.255 


-42) 


8.681 


-42) 


9.117 


-42) 


9.564 


-42) 


1.002 


-41) 



2.405 


'-42) 


2.544 


'-42) 


2.686 


'-42) 


2.833 


'-42) 


2.984 


'-42) 


3.139 


'-42) 


3.299 


'-42) 


3.463 


'-42) 


3.631 


'-42) 


3.804 


'-42) 


3.981 


'-42) 


4.163 


'-42) 


4.349 


'-42) 


4.539 


'-42) 


4.734 


'-42) 


4.933 


'-42) 


5.137 


'-42) 


5.346 


'-42) 


5.887 


'-42) 


6.456 


'-42) 


7.054 


'-42) 


7.681 


'-42) 


8.338 


'-42) 


9.024 


'-42) 


9.740 


'-42) 


1.049 


'-41) 


1.126 


'-41) 


1.207 


'-41) 


1.291 


'-41) 


1.378 


'-41) 


1.468 


'-41) 


1.561 


'-41) 


1.657 


'-41) 


1.757 


'-41) 


1.859 


'-41) 


1.965 


'-41) 


2.074 


'-41) 


2.187 


'-41) 


2.303 


'-41) 


2.422 


'-41) 


2.545 


;-4i) 



1.064 


'-42) 


1.147 


'-42) 


1.233 


'-42) 


1.322 


'-42) 


1.415 


'-42) 


1.510 


'-42) 


1.609 


'-42) 


1.712 


'-42) 


1.817 


'-42) 


1.925 


'-42) 


2.037 


'-42) 


2.152 


'-42) 


2.270 


'-42) 


2.392 


'-42) 


2.516 


'-42) 


2.644 


'-42) 


2.775 


'-42) 


2.909 


'-42) 


3.258 


'-42) 


3.626 


'-42) 


4.015 


'-42) 


4.422 


'-42) 


4.849 


'-42) 


5.295 


'-42) 


5.760 


'-42) 


6.244 


'-42) 


6.747 


'-42) 


7.268 


'-42) 


7.809 


'-42) 


8.367 


'-42) 


8.944 


'-42) 


9.539 


'-42) 


1.015 


'-41) 


1.078 


'-41) 


1.143 


'-41) 


1.210 


'-41) 


1.278 


'-41) 


1.348 


'-41) 


1.420 


'-41) 


1.494 


'-41) 


1.569 


;-41) 
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TABLE I. (continued) Total cross sections for v — d reactions in units of cm^. The 



-X m 



parentheses denotes 10 ^; thus an entry hke 4.279 (-47) stands for 4.279 xlO 



-47 pj^2_ 



E^ [MeV] 



v d ^ V p n 



u d ^ u p n 



Ve d ^ e p p 



Ve. d 



n n 



25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

60 

65 

70 

75 

80 

85 

90 

95 

100 

105 



1.199 


'-41) 


1.317 


'-41) 


1.440 


'-41) 


1.569 


'-41) 


1.704 


'-41) 


1.845 


'-41) 


1.992 


'-41) 


2.145 


'-41) 


2.304 


'-41) 


2.469 


'-41) 


2.640 


'-41) 


2.817 


'-41) 


3.001 


'-41) 


3.190 


'-41) 


3.386 


'-41) 


3.588 


'-41) 


3.796 


'-41) 


4.011 


'-41) 


4.232 


'-41) 


4.459 


'-41) 


4.692 


'-41) 


4.932 


'-41) 


5.178 


'-41) 


5.430 


'-41) 


5.689 


'-41) 


5.954 


'-41) 


6.226 


'-41) 


6.504 


'-41) 


6.788 


'-41) 


7.079 


'-41) 


7.376 


'-41) 


8.957 


'-41) 


1.070 


'-40) 


1.260 


'-40) 


1.465 


'-40) 


1.686 


'-40) 


1.922 


'-40) 


2.172 


'-40) 


2.437 


'-40) 


2.715 


'-40) 


3.007 


;-40) 



1.049 


-41) 


1.145 


-41) 


1.245 


-41) 


1.350 


-41) 


1.458 


-41) 


1.570 


-41) 


1.685 


-41) 


1.805 


-41) 


1.928 


-41) 


2.055 


-41) 


2.186 


-41) 


2.320 


-41) 


2.458 


-41) 


2.600 


-41) 


2.745 


-41) 


2.893 


-41) 


3.045 


-41) 


3.200 


-41) 


3.359 


-41) 


3.521 


-41) 


3.686 


-41) 


3.854 


-41) 


4.026 


-41) 


4.201 


-41) 


4.379 


-41) 


4.559 


-41) 


4.743 


-41) 


4.930 


-41) 


5.120 


-41) 


5.313 


-41) 


5.509 


-41) 


6.528 


-41) 


7.612 


-41) 


8.757 


-41) 


9.959 


-41) 


1.121 


-40) 


1.250 


-40) 


1.383 


-40) 


1.520 


-40) 


1.660 


-40) 


1.803 


-40) 



2.671 


'-41) 


2.933 


'-41) 


3.209 


'-41) 


3.499 


'-41) 


3.803 


'-41) 


4.121 


'-41) 


4.454 


'-41) 


4.802 


'-41) 


5.164 


'-41) 


5.541 


'-41) 


5.934 


'-41) 


6.342 


'-41) 


6.765 


'-41) 


7.204 


'-41) 


7.659 


'-41) 


8.130 


'-41) 


8.617 


'-41) 


9.120 


'-41) 


9.639 


'-41) 


1.018 


'-40) 


1.073 


'-40) 


1.130 


'-40) 


1.188 


'-40) 


1.248 


'-40) 


1.310 


'-40) 


1.374 


'-40) 


1.440 


'-40) 


1.507 


'-40) 


1.575 


'-40) 


1.646 


'-40) 


1.718 


'-40) 


2.107 


'-40) 


2.540 


'-40) 


3.018 


'-40) 


3.540 


'-40) 


4.108 


'-40) 


4.721 


'-40) 


5.378 


'-40) 


6.079 


'-40) 


6.824 


'-40) 


7.612 


;-4o) 



1.646 


'-41) 


1.805 


'-41) 


1.971 


'-41) 


2.143 


'-41) 


2.322 


'-41) 


2.507 


'-41) 


2.698 


'-41) 


2.896 


'-41) 


3.099 


'-41) 


3.309 


'-41) 


3.525 


'-41) 


3.746 


'-41) 


3.973 


'-41) 


4.206 


'-41) 


4.445 


'-41) 


4.689 


'-41) 


4.938 


'-41) 


5.193 


'-41) 


5.453 


'-41) 


5.718 


'-41) 


5.988 


'-41) 


6.264 


'-41) 


6.544 


'-41) 


6.829 


'-41) 


7.119 


'-41) 


7.413 


'-41) 


7.712 


'-41) 


8.016 


'-41) 


8.324 


'-41) 


8.636 


'-41) 


8.953 


'-41) 


1.060 


'-40) 


1.233 


'-40) 


1.415 


'-40) 


1.606 


'-40) 


1.802 


'-40) 


2.004 


'-40) 


2.212 


'-40) 


2.424 


'-40) 


2.640 


'-40) 


2.859 


;-40) 
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TABLE I. (continued) Total cross sections for u — d reactions in units of cm^. The 



-X m 



parentheses denotes 10 ^; thus an entry hke 4.279 (-47) stands for 4.279 xlO 



-47 pj^2_ 



E^ [MeV] 


ly d ^ 


y p n 


ly d ^ 


y p n 


Ue d ^ e p p 


Ve d ^ e^ n n 


110 


3.313 


;-40) 


1.949 


-40) 


8.440 (-40) 


3.081 (-40) 


115 


3.630 


;-40) 


2.097 


-40) 


9.307 (-40) 


3.306 (-40) 


120 


3.958 


;-4o) 


2.247 


-40) 


1.021 (-39) 


3.532 (-40) 


125 


4.298 


;-4o) 


2.397 


-40) 


1.116 (-39) 


3.760 (-40) 


130 


4.648 


;-4o) 


2.549 


-40) 


1.214 (-39) 


3.990 (-40) 


135 


5.009 


;-4o) 


2.702 


-40) 


1.315 (-39) 


4.220 (-40) 


140 


5.378 


;-4o) 


2.855 


-40) 


1.420 (-39) 


4.452 (-40) 


145 


5.756 


;-4o) 


3.009 


-40) 


1.528 (-39) 


4.684 (-40) 


150 


6.143 


;-4o) 


3.163 


-40) 


1.639 (-39) 


4.918 (-40) 


155 


6.539 


;-4o) 


3.318 


-40) 


1.753 (-39) 


5.151 (-40) 


160 


6.941 


;-4o) 


3.472 


-40) 


1.870 (-39) 


5.385 (-40) 


165 


7.350 


;-4o) 


3.627 


-40) 


1.989 (-39) 


5.621 (-40) 


170 


7.765 


;-4o) 


3.781 


-40) 


2.111 (-39) 


5.856 (-40) 



TABLE II. Contributions of the two lowest partial waves. For several representative values of 
the incident neutrino energy E^ are shown the ratios, cr(^5'o)/cj(all) and ^j^q(j(^ Pj) / (j{al\), as 
defined in the text. 



d{v, u)pn 



d{iy, e )pp 



^^[MeV] 
5 

10 
20 
50 
100 
150 



5*0 
0.999 
0.995 
0.972 
0.827 
0.589 
0.433 



0.001 
0.005 
0.027 
0.158 
0.334 
0.410 



^^0 

0.999 
0.993 
0.964 
0.804 
0.561 
0.409 



'Pj 
0.001 
0.007 
0.035 
0.182 
0.366 
0.442 
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TABLE III. Contributions of meson exchange currents to the total cross section. The second 
column (lA) gives the total cross section obtained with the lA terms alone (all the cross sections 
in this table are normalized by the cross sections obtained in our standard run) . The third column 
{+Amex) shows the cross section that includes the contribution of the space component of the 
axial exchange current, while the fourth column {+Akdr,o) gives the results that contain the 
additional contribution of the time component of the axial exchange current. The last column 
{+V'mex) gives results including the full exchange currents using Eq.(^3|) for the vector current, 
i.e., without invoking the Siegert theorem. 



di^v, u)pn 



-^MeV] 


lA 


+Amex 


5 


0.949 


1.000 


10 


0.942 


0.999 


20 


0.934 


0.996 


50 


0.927 


0.991 


100 


0.925 


0.984 


150 


0.924 


0.979 



A^KDRfl 


+y'MEX 


0.999 


1.000 


0.999 


1.000 


0.996 


1.000 


0.991 


0.999 


0.984 


0.997 


0.979 


0.996 



d{v, e )pp 



-4MeV] 


lA 


+AmEx 


+Akdr,o 


5 


0.952 


0.999 


0.999 


10 


0.945 


0.997 


0.997 


20 


0.937 


0.994 


0.994 


50 


0.928 


0.985 


0.985 


100 


0.924 


0.974 


0.974 


150 


0.922 


0.966 


0.966 



+n 



MEX 

1.000 
1.000 
1.000 
0.999 
0.995 
0.993 
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TABLE IV. Model dependence of total cross sections. The second column (YHH) and the 
third column (KN) give a{Y}i}i)/ a {standard run) and a(KN)/ a {standard run), respectively. 
The fourth column (Reid93) [fifth column (NIJ II)] gives the ratio of the total cross section ob- 
tained with the use of the Reid 93 potential [Nijmegen II potential] to that of our standard run. 
The last column (A(CRSW)) gives the ratio of the total cross section obtained with the A-current 
of Carlson et al. to that of our standard run. 



d{u, u)pn 


E^[Me\] 


YHH 


KN 


Reid93 


Nun 


A(CRSW) 


5 


0.962 


1.002 


0.997 


1.002 


0.965 


10 


0.955 


1.003 


0.998 


1.002 


0.961 


20 


0.946 


1.000 


0.998 


1.001 


0.956 


50 


0.964 


0.993 


0.999 


1.000 


0.953 


100 


0.961 


0.971 


1.000 


1.000 


0.953 


150 


0.915 


0.943 


1.000 


0.999 


0.954 


d{iy,e~)pp 


K[MeV] 


YHH 


KN 


Reid93 


NIJ II 


A(CRSW) 


5 


0.956 


1.019 


1.003 


1.003 


0.968 


10 


0.949 


1.008 


1.003 


1.002 


0.964 


20 


0.948 


1.002 


1.002 


1.001 


0.959 


50 


0.961 


0.990 


1.001 


1.000 


0.956 


100 


0.955 


0.968 


1.001 


0.999 


0.956 


150 


0.897 


0.941 


1.001 


0.999 


0.956 
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TABLE V. Model dependence of i? = a{NC)/a{CC) = a{vd -^ vnp)/a{ued -^ e'pp). For 
representative values oi Ey, R for our standard run is given in the second column. The third 
through the sixth columns give Ra^ with a = lA, Reid93, NIJ II, and A(CRSW), normalized by 
Rstandardrun- See also the Caption for TABLE ITV. 



;[MeV] 


^standard run 


lA 


Reid93 


NIJ II 


A(CRSW) 


5 


0.277 


0.997 


0.994 


0.999 


0.997 


10 


0.410 


0.997 


0.996 


1.000 


0.997 


20 


0.447 


0.997 


0.997 


1.000 


0.997 


50 


0.433 


0.999 


0.998 


1.000 


0.997 


100 


0.398 


1.001 


0.999 


1.000 


0.997 


150 


0.375 


1.003 


1.000 


1.001 


0.998 
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FIGURES 
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FIG. 1. Total cross section for radiative neutron capture. The solid curve corresponds to the 
results of our full calculation including the lA and exchange currents and all the multipole ampli- 
tudes. The dashed and dash-dotted curves show the individual contributions of the niagnetic-dipole 
and electric-dipole amplitudes, respectively. The data are taken either from the neutron capture 
reaction itself |^^, or from its inverse process |48,49|, with the use of detailed balance for the latter. 
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FIG. 2. Total cross section for the reactions: Ved -^ e~pp, Ved -^ e~^nn, I'd -^ upn, and 
i?d -^ Upn. The sohd and dotted curves show the charged-current reaction cross sections for v and 
i>, respectively, while the long-dash and dash-dotted curves give the neutral-current reaction cross 
sections for v and z/, respectively. 
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FIG. 3. Electron energy spectra for the i^ed -^ e pp reaction. 
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FIG. 4. Electron angular distribution for the fed -^ e pp reaction. The solid curves show 

the results of our standard run, while the dotted curves correspond to the simplified expression, 

Eq.([69|), normalized to the standard run results at 9l = 0. 
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FIG. 5. Neutron energy spectra for the lyd -^ i/pn reaction. 
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FIG. 6. Neutron angular distribution for tlie vd -^ vpn reaction. 
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FIG. 7. Kinematically allowed region of the electron scattering angle 6l in the U(,d -^ e~pp 

reaction at E^ = 150 MeV. The dotted area represents the allowed region. The constraint on 6l 

sets in at EL = E'f. 
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FIG. 8. Kinematically allowed region of p', the relative momentum of the final two nucleons in 



the t'ed -^ e pp reaction at Eiy = 150 MeV. The lower limit of p' reaches zero at E'^ = E[ 
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